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Abstract: In this paper we present a self consistent theoretical approach for the calculation of the Survival 
Probability for central dijet production . These calculations are performed in a model of high energy soft 
interactions based on two ingredients: (i) the results of N=4 SYM, which at the moment is the only theory 
^ that is able to deal with a large coupling constant; and (ii) the required matching with high energy QCD. 
C Assuming, in accordance with these prerequisites, that soft Pomeron intercept is rather large and the slope 
! of the Pomeron trajectory is equal to zero, we derive analytical formulae that sum both enhanced and 
^ semi-enhanced diagrams for elastic and diffractive amplitudes. Using parameters obtained from a fit to the 
K> available experimental data, we calculate the Survival Probability for central dijet production at energies 
^ accessible at the LHC. The results presented here which include the contribution of semi-enhanced and net 
. 5^ (diagrams, are considerably larger than our previous estimates. 
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1. Introduction 

The original evaluation of the survival probability for a large rapidity gap (SPLRG) was made by Bjorken 
[l](see also Ref. [2]), in a one channel formalism, where the elastic amplitude was assumed to have a 
Gaussian shape in impact parameter space. This enabled him to obtain a closed analytic expression for 

< |S 2 | >, the SPLRG. 

We [3] extended Bjorken's treatment to five different models, all of which reproduce the measured total 
cross sections at energies upto and including W = 1.8 TeV. However, the values obtain from these models 
for < \S 2 \ > differed greatly, e.g. for the Tevatron energy they predicted values between 9.6 and 32.6%. 
Indicating that using parameters determined only from the elastic amplitude, was insufficient to determine 

< |5 2 | > uniquely. 

We then broadened our approach to a two channel Good- Walker description, which included elastic as 
well as diffractive channels [4]. This was called for, as the experimental data from the Tevatron [5] indicated 
that the diffractive channels were important and should not be neglected, asd+add m 0.85. Including the 
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diffractive channels necessitated including also the triple Pomeron coupling, which plays a crucial role in 
determining < \S 2 \ >. 

The difficulty inherent in including the diffractive channels (i.e. 2>JP coupling), is that one has to sum 
over all orders (loops) of these interactions, and there is no rigorous method for evaluating the so called 
enhanced and semi-enhanced (Pomeron) diagrams. 

To make the problem tractable, we have in two previous papers (see Refs. [6,7]) constructed a model 
of soft interactions at high energies based on the postulates that stem from N=4 SYM and QCD. The 
model includes the following features: 

1. It is built using Pomeron and Reggeons as the main ingredients as follows from N=4 SYM( see 
Refs. [8-12]); 

2. The intercept of the Pomeron should be rather large. In N=4 SYM [8, 13] we expect Ajp = 
ajp(O) — 1 = 1 — 2/y/X ~ 0.11 ~ 0.33, based on estimates of A = 5 ~ 9 derived from the early LHC 
measurements from the cross section of multiparticle production as well as from DIS at HERA [14]. 

3. As follows from N=4 SYM, the slope of the Pomeron «jp(0) = 0, this is compatible with our fits 
to the cross section data in the ISR-Tevatron energy range [6,7]. 

4. A large contribution to the cross section comes from processes of diffraction dissociation, as in N=4 
SYM at large A only these processes contribute to the scattering amplitude. In other words, in this 
model the Good- Walker mechanism [15] is the main source of diffractive production; 

5. The Pomeron self-interaction should be small (of the order of 2j\f\ in N=4 SYM), and much 
smaller than the vertex of interaction of the Pomeron with a hadron, which is of the order of A; 

6. The last condition is not a prerequist of N=4 SYM, but follows from the requirement of the natural 
matching with perturbative QCD, where the only vertex that contributes is the triple Pomeron vertex 
[16]. This ingredient differentiates our model from other models on the market (see Refs. [17, 18]); 

In this paper we apply the theoretical approach that was developed to calculate the value for SPLRG 
processes at the LHC. This problem has been a subject of much discussion [6,17] since our estimates first 
given in Ref. [6], suggested a small value of the survival probability for central Higgs boson production, 
in disagreement with the estimates of the Durham group [17]. In Ref. [6] we only included the enhanced 
Pomeron diagrams. In this paper we present our estimates for the full set of the diagrams. We discuss 
the important topic : the derivation of formulae for SPLRG processes originating from three sources: the 
Good- Walker mechanism, semi-enhanced diagrams and enhanced diagrams for the Pomeron interaction. 
We discuss this derivation in the first three sections, while the fourth section is devoted to numerical 
estimates. 
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Figure 1: Survival probability for exclusive central diffractive production o_f the Higgs boson. Fig. HTa shows the 

| &P~ Q 1 . 3. 

contribution to the survival probability in the G-W mechanism, while Fig. |l|-b illustrates the origin of the additional 
factor (| Sl nh |). 

2. The main formulae: Good- Walker mechanism 

The calculation of the survival probability for the hard processes is a choice example that shows that 
the cross section of the hard processes cannot be calculated, except for some special cases (total DIS 
cross section, inclusive cross sections), without a substantial knowledge of the 'soft' physics. The general 
formula for the Good- Walker mechanism (sum of the non-enhanced Pomeron diagrams) and for the sum 
of enhanced diagrams has the form 

/ dH fc fc <p^ > 2 <p\k > 2 e -K*W) Sfj k (b)} 2 

(\S 2 \} = l — ' 2, (2.1) [sp 

/ d2b {Ei.fc <P\i > 2 <p\k > 2 / d 2 b'A i H {s,b)A k H (s,b-b')j 

where, < p\i > is equal to {^proton \ ^i) and , therefore, < p\l >= a and < p\2 >= (3, a 2 + fi 2 = 1. Our 
notation as well as the derivation of Eq. ( pTl| ) can be found in Ref. [19]. 

The factors exp[— are responsible for the Good- Walker mechanism contribution for the survival 
probability, while the factors S^ k (b) account for an additional suppression due to a structure of the Q% k- In 
Fig. Ifwe illustrate examples of the different contributions in S^ k . Fig. jj)^ which is taken from our paper [6] 
shows that S^ k depends critically on both the hard and the soft amplitude. Our main goal is to find this 
dependence and calculate S^ k . The parametrization of the hard amplitude A l H (s, b{) has been discussed in 
Ref. [6], and we refer the reader to this paper for details. 

3. Enhanced diagrams 

The calculation of the sum of the enhanced diagrams has been treated in detail in our previous papers [6,7] 
and so we will not repeat the discussion here. 
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di-jets, Higgs 
with mass M 



b) 





G-W contribution 



c) 




G-W contribution 



Figure 2: The set of diagrams that is selected and summed for the calculation of the survival probability for 
diffractive Higgs production. Fig. pfa shows the diagrams in G-W + enhanced diagrams approach, in Fig. p^b the 
same approach is shown, however, to calculate the value of the survival probability we add the first semi-enhanced 
diagram. The approach for giT(Y) w 1 but ApT(Y) <C 1 (net diagrams) is shown in Fig. pfc. 
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In the case of summing the enhanced diagrams the damping factors Sf k {b; 6162) have a simple structure 



SP k (b;hb 2 ) = (\S enh (Y)\) A i H (s,b 1 )A k H (s,b 2 ), 



(3.1) 



SPED1 



where (| S en h (Y) |) depends only on energy (total rapidity Y = ln(s/so)) and it can be considered to be a 
common factor. Although, this factor appears in Ref. [6] , for completeness of the presentation it is given 
below. 
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(\S 2 enh (Y)\) 



where 



S(T) 



S(T(Y)) - 2e~ A ^ Y ~ yh ^ 2 S1(T(Y)) + e -^(X-vh)/2 S2 ; 



Sl(T) = ^l-T{l + T) + (l + 2T)eTr(o,i 
S2(T) = ^ |t[(T-1) 2 -2] + (l + 3T) e *r(().- 



r(y) = 7 ( e A -( y - y ') - l) (e A - r ' - l' 

r (0, l/T) is the incomplete gamma function (see formulae 8.35 in Ref. [20]). 

4. Enhanced diagrams and semi-enhanced diagrams as a perturbation 



(3.2) 

(3.3) 
(3.4) 
(3.5) 

(3.6) 



TT 







a) 



b) 



Figure 3: The set of diagrams in terms of the exact Pomeron and exact vertices that we include in our perturbative 

approach for the summation of semi-enhanced diagrams. 

(seen „ 
Fig. g-b shows the changes in S" k . 



B een 
-a gives contribution to fi^fc (see section 2.3) while 



As has been discussed in our approach based on the summation of the enhanced diagrams (see Ref. [6]) 
the contribution of the semi-enhanced diagrams to the value of the survival probability has been neglected. 
The first correction due to such a contribution is shown in Fig. j^b ( compare with Fig. j^a). The expression 
for this correction has the form (see Fig. ^-bj 



(\S£ k (Fig.^pb;Y,b)\) = J d 2 h { (| S 2 enh (Y) |) Af (Y, b x ) Af (Y, b-b x 
rvi ryi 

2 1 d Vl dy 2 / dy 3 G (Y - yi ) T ( Vl ,y 2 , y 3 ) G (y 2 - 0) <| S 2 enh (y 2 ) |) 



x A" [Y,b- b x 



d 2 b'Af (Y,b') ~g k [bi 



(4.1) 



(4.2) 



The coefficient 2 stems from two diagrams where the exact Pomeron is attached to upper and lower proton 
in Fig. H^bT G(Y) is the exact Green's function of the Pomeron and is given by 



G(Y) 



1 



exp 



T (Y) J T(Y) 



r 



T(Y) 



(4.3) 
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with 



T(Y) = 7 e" 

r (2/1)2/2)2/3) is the exact vertex that has been discussed in section 3.4 of Ref. [7]. It is equal to 



r (s/1,2/2, J/3) 



1 



, , r , r{r 1 (2T(y 1 -2/ 2 )) - r 1 (2T(y 1 -2/ 3 ))} 

4 T(yi - y 2 ) - T(y 1 - y 3 ) 

withTi(T) = (1/T 3 ) x{T(l + T)-exp(-l/T) (1 + 2T) r(0,l/T)}. 

To calculate the survival probability we need to introduce the additional change in 0^ ^ in Eq. 
(see for example Ref. [21] for detailed arguments). The resulting & has the form (see Fig. ^-&) 



(4.4) | ES11 

(4.5) 
(4.6) 
(4.7) 



fyfc (Y,b) 



yi 



d 2 b 1 [G(Y)~g i (61)^(6-^) 

/Vl 
dy 3 G(Y - y x ) V { yi ,y 2 , y 3 ) G {y 2 - 0) G (y 3 - 0) 



5. Net diagrams. 
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(4.9) 




TV 



G(F,y) 



a) 



pad 



G(F,y') 

6) 



G(Y',2/) 



Figure 4: The set of diagrams( see Fig. ^-a) for S^, in the approximation where gY (Y) ss 1 while AT (Y) <C 1 ( 
net diagrams) and the graphic form of the equation for S^ k = G^k- Zigzag line denotes the hard Pomerons. 

In this section we discuss the survival probability in the approximation where g(F (Y) ~ 1 while 
AjpT (Y) <C 1. In this approximation we need to sum a 'net' system of the diagrams, an example of which 
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is shown in Fig. |pa7 To find S^ k we introduce the function G^ (Y; y) for which we can write the equation 
(see Fig. flflj 

G(Y,y) = A%(h) - A J V dy'G{Y,y') G (Y, y') . (5.1) 

Eq. ( |b.l[ ) can be rewritten in the differential form as 

dG{Y,y) 



dy 



AG{Y,y') G (Y, y') 



(5.2) 



eqnd 



SPND1 



SPND2 
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The solution to Eq. ( p. 4 ) is 

G(Y,y) = Ag(h) exp 



AG(Y,y') dy' 



1 



1 + gi(b 2 )T(Y) + ~g k (h) T (y) 



(5.3) 



S^ k can be expressed with the aid of the functions G in the following way: 



S H (Y,b) = [fhAftf-ftAgMx- T-^-- " 

J (l + (G 3JP /7) [9k(b ~ h) T (Y) + ~ gi (bi) T (Y/2 - y h /2) 



1 

l + {G ZIP h)(Ub-bi)T{Y) + ~g k (h)T(Y/2 + y h /2) 



(5.4) 



where yh = ln(M 2 / sq) and M denotes the mass of the produced system (dijets, Higgs boson etc.) and 
s = 1 GeV 2 . 

Finally to calculate the survival probability we require Qi tk which was derived in Ref. [7] 







1 



i k 



■2iJ 



d Z b 



9i(b') ~g k (b-b>) (I/7) T(Y) 



1 + (Gsip/7) T(Y) 



~gi(b')+9k[b-b' 



(5.5) 



SPND5 



6. Net and enhanced diagrams. 





jjqnden 

Figure 5: The net and enhanced diagrams (Fig. H-a) and the set of diagrams with the exact Green's function of 

'"qnden 



Pomeron and exact vertices (Fig. y-b). Zigzag line denotes the hard amplitudes. 

In the spirit of our approach we can generalized Eq. ( p™^ and Eq. ( Ib^^ taking into account the 
enhanced diagrams. As we have discussed, we need to replace T (Y) by G (T (Y)) in these equations (see 
Fig. |hPa and Fig. |^bj and to multiply Eq. (p^J by (| S^ nh (Y) |). This factor is responsible for summing 



eqnden 



the diagrams where the Pomeron lines are attached to the hard amplitude ('hard Pomeron') (see Fig. Epaj . 
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7. Numerical estimates and prediction for the LHC 



W=1.8TeV 

1.96 TeV 
W=2.36 TeV 




M 2 (MeV 2 ) 

Figure 6: The estimates for the value of SP for dijets with mass M. The dotted line corresponds to the Higgs boson 
production with Mhi ggs = 100 GeV . We estimate a 60 4- 70% margin of error in our results. 

Using Eq. ( p]i| j and Eq. (j§7^T~we have calculated the survival probability as a function of produced 
di-jets mass at cm. energies of interest. We estimate the margin of errors of our results to range from 60% 
at W=2 TeV to 70% at W=14 TeV. The central values obtained are shown in Fig. 

The main feature of our SP results is that their values are considerably larger than the values obtained 
in our previous studies [7]. Our current estimates are compatible with the CDF Tevatron data [22] as well 
as with the theoretical estimates [23,24]. 

Our calculations presented in Ref. [6] only summed over the enhanced diagrams, while our present 
calculations sum over the complete set of enhanced, semi-enhanced and net diagrams. Comparing the 
fitted values of our parameters in the two studies we note that our present fitted values of 7 and G^p 
are much smaller than those obtained in Ref. [6]. Consequently, the screening initiated by the enhanced 
diagrams with the present parameters is much smaller that the screening caused by the same diagrams in 
our previous study. On the other hand, in our both papers we found out that (i) the main source of the 
small values of the SP is the enhanced diagrams; and (ii) the enhanced diagrams are the only cause for the 
dependence of the SP value on the mass of produced di-jets. Therefore, the rather large value of SP and 
its mild dependence of the mass of produced di-jets have the same origin: the small values of G^jp. 
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8. Conclusions 

Fig. |§f shows our estimates and predictions for the value of the SP. These values are the results of our 
theoretically self consistent approach that has been discussed in Refs. [6,7]. We predict the value of the 
SP for Higgs production {Mniggs = 100 GeV) of about 3 — 5% for the LHC range of energies. These values 
are much larger our previous estimates, for a model in which only enhanced diagrams were taken into 
account [6]. Our main conclusions are that the value of the SP, as well as its mass dependence are very 
sensitive to both the particular form adopted for the Pomeron interaction, and to the values of the fitted 
parameters that determine the strength of the contribution of the enhanced diagrams. 

The lesson from our present study is straight forward and discouraging to an extent, since we need a 
thorough knowledge of the theory of interacting Pomerons to guarantee the value of the SP. On the other 
hand, the model we constructed incorporates the main features of two compelling theoretical approaches: 
N=4 SYM and perturbative QCD. 

We conclude that our theoretical knowledge of Pomeron interactions, as well as the set of available 
experimental data, are not sufficient to determine the strength of the Pomeron interaction with adequate 
accuracy, so as to provide precise estimates for the SP. 
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